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I have been asked to provide an elementary introduction to phase transitions and 
critical phenomena. The subject is vast; thus, this can really only be a somewhat superficial 
personal overview of the subject. Many important topics are left out; at the end of these 
notes is a brief bibliography of sources for further reading. 

While most of the universe is a dilute gas of hydrogen, phase transitions are crucial 
to our existence. We breath air, blood flows through our veins, and our bones are solid. 
Water boils and lakes freeze over. Mathematically, however, phase transitions are rather 
remarkable. Statistical mechanics is based on Boltzmann weights 

e -f3E(S) 

P(S) = — 

where the inverse temperature /? = pf • From these the partition function 

0E(S) 



e 



is constructed. But e~@ E is an analytic function of 0, i.e. it has no singularities. 

Phase transitions can only occur when an infinite number of states are available, such 
as with infinite volume V. Start with the partition function 



E 

where N(E) is the number of states of a given energy E. Each piece of a large system can 
have its own energy, so one should expect (E) proportional to V. Work with the energy 
density £ = E/V. Since there are of order V places to put each bit of energy, we also 
expect the number of states of a given energy density to grow exponentially 



N(E) = exp(VS(S)) 
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This defines the entropy density S at the given energy density. Pulling out the volume 
factors explicitly, 

Z = V J d£e v W-M 
As V — > oo the integral is dominated by the maximum of the integrand, where 

or the more usual form 

AS = AQ/T 

In this saddle point approximation 

Z ~ e~ t3VF 

where free energy is 

F = S — S/(3 

and all of thermodynamics follows. 

A singularity at a phase transition requires structure in N(E). Such structure also 
requires spatial correlations. Otherwise 

z(v) ~ z(v/n) n 

and the free energy becomes just that for a small volume. 

In much of the following I will use the Ising model as an example. This has "spins" 
Si E {1, —1} occupying lattice sites and has energy E = — s i s ji where {ij} denotes 

nearest neighbors. Table 1 gives N(E) for the two dimensional Ising model on 10 by 10 
lattice with periodic x boundaries, cold y walls. The length of the line gives the entropy. 
The phase transition is hidden in a subtle flatness. 

The numbers in this table are rather large. They add to 2 100 = 1.267 x 10 30 . Our 
universe is only about 10 27 nanoseconds old, suggesting that it is impossible to calculate 
this list by simple counting. This is a frequent argument for Monte Carlo. See my state- 
counting papers to learn how I got this list. 

How can we get enough correlation for a phase transition? One way is surface tension. 
Let small volumes be in two possible phases, i.e. water and steam. Suppose we pay a 
penalty for an interface between the phases. For a model, put the system on a lattice with 
two states in each cell, Sj G {1,-1}, representing water and steam 

Z = J2e~^ F(s ' )+ ^ Js ' S3 . 

Si 

Here J represents the surface tension and nn means nearest neighbors. The relative free 
energies of the two states can be influenced by, say, the pressure. The interesting case is 
when they are near each other, so let me expand 

F(s i ) = F-Hs i 
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Table 1. State counts for the Ising model on a 10 by 10 lattice. 
ff/2 + 100 N(E) E/2 + 100 N(E) 
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Figure 1. Phase diagram for the Ising model. 



Figure 2. Simulation of boiling water. 

where F = (F(l) + F(-l))/2 and H = -(F(l) - F {-!))/ 2. 

Thus we expect the boiling of water to be described by something like the Ising model 
in an applied field, more usually used for magnets 

Z = J2e^ n ^^+ H E i ^ 

Si 

Boiling occurs on a flip in the sign of H at low temperature. The basic structure of this 
model, shown in Fig. 1, is a first order line running from a critical f3 to infinity along the 
(3 axis. The H —H symmetry, expected for magnets, is broken in water by higher order 
effects. 
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Fig. 2 shows a simulation of water boiling, obtained by taking an ordered Ising system 
and turning on a field in the opposite direction to the magnetization. The bubbles nucleate 
from larger fluctuations. The picture was made using xpotts, a Potts model simulator from 
my xtoys collection at [http: / / penguin.phy.bnl.gov / www / xtoys / xtoys. html . 



OK, so we need infinite volume, but how infinite? Try a chain of sites in a ring, let 
the ring length go to infinity. A bond in the chain connecting s and ^'contributes 



- \ ( -0 e P-H 



T , = e P*s'+H(s+s')/2 



s,s' 



This is the "transfer matrix." Summing over spins gives the partition function for an N 
site lattice 

Z = Tr(T N ) 
This can be calculated by diagonalizing T 

Z = A? + A^ 

with 

A± = e p (cosh(H) ± ^/sinh 2 (ff) + e" 4 / 3 

The free energy per site is dominated by the largest eigenvalue 

-0F = ^ = log(A+) + exp(-iVlog(A + /A_)) + . . . 

Since A_|_ is analytic and positive, the infinite volume free energy has no singularities. We 
fail to find a phase transition. Note that the finite volume corrections are exponentially 
small. This shows that the theory has a mass gap. As (3 goes to infinity with H = 0, 
the eigenvalues become equal and the mass gap goes to zero. The phase transition, to the 
extent there is one, occurs at zero temperature. 

Physically, a kink anti-kink pair has a finite probability, but an infinite number of 
possible separations. This infinity always disorders the system. In more dimensions a 
big bubble pays a penalty proportional to its surface, so it is suppressed. For now I will 
continue to concentrate on the H = Ising model, but with more neighbors. 

So we need more neighbors. Now for an amusing model. Suppose each site has 2d 
neighbors, but on a system without closed loops. Let me build this lattice recursively. A 
single outermost site is connected to a deeper neighbor, which has 2d-l other neighbors. 
This site can have either spin -1 or spin 1. Fix this spin and define Z(s) to be the partition 
function obtained by summing over all deeper spins. This partial tree is then recursively 
built up, defining a Cayley tree as sketched in Fig. 3. To get things started, apply an 
infinitesimal field to the zeroth level sites 



Z (s) = e 
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Figure 3. Constructing a Cayley tree. 




1 2 3 4 5 

Figure 4. Recursive solution of the Ising model on a Cayley tree. 



At level n, the partition function can be written as a sum over the values of the level n — 1 
neighbors 

Z n (l) = (e^ n _ 1 (l) + e-^ n _ 1 (-l)) 2d - 1 
Z B (-1) = [ePZ n _ x {-\) +e-?Z n _ 1 (l)) 2d - 1 



Dividing we obtain 



Rr 



ZJl) 



n-1 



Z n (-1) \eP + e~PR 



ln-1 



2d-l 



We can now look for a stable asymptotic solution for a nontrivial R by asking that 
R n = Rn-i = R- As shown in Fig. 4, R = 1 is always one solution, but more appear at 
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Figure 5. A two dimensional lattice and its dual. 

the critical coupling which occurs when the derivative of the right hand side is unity. This 
happens at 

tanh(A:) = ^-j 

When (3 exceeds the critical value, a value of R above one will flow towards the non-trivial 
fixed point as seen in the figure. The transition for this model is second order since the 
fixed point smoothly moves to unity as the critical coupling is approached. This concept 
of iteration generating a "flow" towards a "fixed point" will be a recurring theme in the 
following. Note that as d goes to one the critical temperature moves to infinity. One 
dimension is a critical case. 



I now digress into the topic of duality, which gives the exact critical point for the 2d 
case. To see how this works, change variables from sites to bonds. This can be done a 
couple of different ways. 

For a given configuration, each bond is either excited or not. An excited bond gives 
a factor of e - ^. A non-excited bond gives '. Re-express the sum over sites as a sum 
over bonds being excited or not. These variables over-determine the spins, and thus are 
constrained; in going around a loop one must encounter an even number of excited bonds. 
The product of bonds around a plaquette is positive. This view of the system forms the 
basis for the low temperature expansion, i.e. at low temperature most bonds not excited. 

Now for an alternative view, write e l3siSj = cosh(/3) + SiSj sinh(/3). On each bond 
either take one or the other of these two terms. Assign 1 to bonds with the first choice, — 1 
to the others. Given a particular configuration of such choices, sum out the spins. This 
will give zero unless there are an even number of bonds coming out of a site which use the 
sinh(/3) term. This expansion on the bonds is the basis of the high temperature expansion, 
i.e. small f3 means small sinh(/3). Again we get a sum over two states for each bond, with 
a constraint on these choices. The bond variables multiplied over the four links coming 
from any bond must be positive. 

These two constrained systems are closely related. Transfer the variables from the 
bonds to dual ones, new bonds crossing the middle of each old one, new sites in the middle 
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Figure 6. Recursively solving the mean field equation. 

of the old plaquettes, as shown in Fig. 5. The product of old bonds out of a site becomes 
a product of new bonds about the new plaquette. With this mapping on one of the bond 
descriptions, the two constrained systems now satisfy the same constraint. This means 
that up to irrelevant factors, the physics at two values of (3 is related 

e~ 2f3 = tanh(/3') 

Under this mapping strong and weak limits interchange. If there is only one singularity, 
it must occur at the self dual point (3 = (3' = (3 C = \ log(l + V2) = .44068 . . .. This 
idea of duality extends to Z n clock models, Potts models, 4-d Z<i lattice gauge theory, 3d 
gauge-Higgs system, QED with monopoles, . . . Seiberg .... But it all starts here. 

I now turn to mean field theory. Consider a large number of dimensions. Then each 
site has lots of neighbors, suggesting we might assume their effect can be averaged. Suppose 
we are in a magnetized state with (s) = M. Look at one spin in the average field of the 
others. Calculate its magnetization 

e 2dMf3 _ e -2dMf3 

M = e 2 dM p + e -2 dM p = tanh(2rfM/3) 
This can be solved by iterating 

M -> M + AM 

AM = tanh(2dM/3) - M 

We have a non-trivial fixed point only if 





ICtl II 1 \.\J A) A 

tanh(x)-x 
tanh(2*x)-x 
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Figure 7. The effective potential from Jensen's inequality. 



As shown in Fig. 6, this "flow" like picture gives a similar result to the Cayley tree case at 
large d, but the approximation misses the transition moving to infinity as the dimension 
goes to one. 

Jensen's inequality provides a more formal approach to mean field theory. Note that 

,2 

e x is a convex function, -^e x = e x > 0. If a; is some stochastic variable, this means that 
(e x ) > e^ x ' . Show this for homework. 

Consider a "fake" weighting with P(si) = e H +e - H ■ With this probability distribution 



Thus H might be thought of as a "source" pulling on the spins. Now I manipulate the 
partition function 



(si) p = tanh(if) 






> exp(((3^2siSj -H^ Sl + yiog(2cosh(iT)))p) 



{ij} i 
= exp(V(df3 tank 2 (H) - H tanh(#) + log(2 cosh(ff)))) 



Thus for any H the true free energy is less than 



F <= F mf = -d[3 tank 2 (H) +Htanh(H) - log(2cosh(if)). 
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Figure 8. Mean field theory for the three state Potts model. 

This is an "effective potential" which, depending on the value of (3, can have one 
minimum at H = or a double well shape with two minima, as sketched in Fig. 7. The 
latter represents the ordered phase. The critical point occurs when 0(H 2 ) term vanishes, 
i.e. the "mass term." This happens at 

d(3 c -1/2 = 

or (3 C = ^j, as before. 

These transitions have been second order, i.e. the system evolves continuously with the 
coupling parameters. The picture is a bit different with three states, where cubic terms can 
drive us to first order transitions, and physics becomes discontinuous. For example consider 
the three state Potts model, a system with a Z3 symmetry. Take S{ £ {1, e 27 ™/ 3 , e -27 ™/ 3 }. 
Suppose the bonds have low energy if the spins are "parallel" or equal, higher energy 
otherwise 

E = -"£Res* Sj 

ij 

Define the magnetization to vanish for a random distribution 

M = (Res,) 

For mean field theory, we are to replace neighbors with the average and solve self consis- 
tently. Now there are two anti-parallel cases, each giving minus half a unit of magnetization 

e d/3M _ e -dpM/2 
M = 

e df3M _|_ 2 e ~df3M/2 

Expand the right hand side for small M 

M = ^/3M + |/3 2 M 2 + 0(M 3 ) 
2 8 
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For the Ising case there was no 0(M 2 ) piece, and the 0(M 3 ) piece was negative. As 
shown in Fig. 8, now a non-trivial solution appears before reaching (3 = 2/d. The new 
solution appears discontinuously; when it becomes of lower energy, we have a first order 
transition! Note that when the extra solution first appears, it is meta-stable and not the 
lowest energy; one can use the Jensen inequality argument to estimate when it drops below 
the unmagnetized state. Simulation results show that for three and more dimensions the 
transition is indeed first order. In two dimensions, however, it is second order and the 
prediction fails. In one dimension there is no transition, just as for the Ising case. 

Lattice gauge theory is summarized by the path integral 

z = e MP ReTrU p) 

U P 

with Up = U1U2U3U4, and the Ui are link variables running around the plaquette in 
question. The local symmetry Uij — ■> giUijgj implies there is no barrier to twirling a local 
group of links around. Without gauge fixing U cannot have an expectation value for any 
p. This is Elitzur's theorem. One should play the Jensen game for more rigor, but proceed 
naively anyway, trying to find a self consistent expectation for a link. Do Z 2 for simplicity, 
which gives the same result as the above Ising case except for the replacement M — > M 3 
for the average field 

M = tanh(/3M 3 ) 

Now there is no linear term at all on the right hand side. The prediction is for a strong 
first order transition. Most lattice gauge transitions in fact are first order: Z2-4 i n four 
dimensions; all known gauge groups in 5 or more dimensions. However, in 3-d, the Z2 
gauge model is dual to the Ising model; so, the transition is second order. In 2-d, gauge 
fixing turns a gauge model into a one dimensional spin system, with no transitions for any 
finite dimensional group. 

Later in these lectures I will generalize these arguments to suggest a first order de- 
confining transitions for pure SU(3) gauge theory at finite temperature and for the chiral 
transition with three massless quarks. This will involve some mathematical formalism that 
I postpone. 

Now I change the subject a bit and remind you of the formal connection between path 
integrals and statistical mechanics. This is one reason quantum field theorists are also 
interested in phase transitions. Let me start with a simple quantum mechanics problem 
with the Hamiltonian 

H = p 2 /2 + V(x) 
Here p and q are conjugate variables with [p, q] = —i. Look at 

Z = Tie-? 11 

i.e. zero space dimensional quantum statistical mechanics. As (3 — > 00 we project out the 
ground state and get ordinary quantum mechanics at zero temperature. This is also the 
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trace of the evolution operator e~ ltH for imaginary time t = —ij3. Quantum mechanics in 
imaginary periodic time is the same problem as quantum statistical mechanics. 
I now break up (3 into a large number of "imaginary time slices" 



N 

-/3H/N 



e 

i 



N 



Insert a complete set of states at each slice 

Z = I d Xl ...dx N flixile-^^lxi + 1) 

i=l 

where xn + i = x\. Now for large iV we can approximate 

( Xi \e-^ H / N \x i+1 ) ~ e-P y ^/ N ( Xl \e- Pp2/{2N) \x l+1 ) 
The second factor can be worked out by inserting a complete set of momentum eigenstates 

( Xi | e -/3p 2 /(2A0| Xj+l) = J dpe ^ P (^ +1 ) e -P P V(2N) = y| e -7V(x l+1 -x I ) 2 /(2/3) 

Thus we have the simple form 

Z = J dxi . . . dxN£~ S 
where the "lattice action" is simply 



and the lattice spacing a = (3/N. This defines the path integral, which formally is a classical 
statistical mechanics problem in one dimension. It represents the thermal dynamics of a 
polymer, with the Xi denoting the coordinates of the atoms in the chain. 

Going from a single particle to a field theory, D space dimensional quantum mechan- 
ics is equivalent to d = D + 1 dimensional classical statistical mechanics. The infinite 
"time" limit gives the ground state, while finite imaginary time relates quantum statistical 
mechanics to classical statistical mechanics in one more dimension. 

Second order phase transitions are essential to continuum limits. Taking the lattice 
spacing to zero requires physical correlation lengths to diverge in lattice units. The particle 
physicist's e~ mr corresponds to the statistical mechanic's e~ n ^. With lattice spacing a, 
we identify r = na and £ = m/a. A continuum limit requires £ — > oo, which occurs at a 
phase transition, most particularly at a second order one. For lattice gauge theory, 4-d is 
a borderline case and the transition occurs at (3 ~ 1/g 2 = oo, with the approach given by 
standard asymptotic freedom arguments. 
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Now for another jump in subject. Despite the usual pedagogical approaches, effective 
potentials want to be convex. In field theory language, consider 



Z = J #e" s ^ 

Adding in some external sources 

Z( J) = J dcf>e- s W +J <t> 

general correlation functions can be found by differentiating with respect to J. Here I 
shorthand notate J<p = J dxJ(x)<p(x) in the continuum, or J<p = Jj0j on the lattice. 
Think of J as an external force pulling on the field. Such a force will give the field an 
expectation value 

where I define the free energy F(J) = — log(Z(J)). Now imagine inverting this to find 
what force </(<&) gives some desired expectation value, i.e. solve 

*(J) = <*>J<») = -% 
In terms of this formal solution, construct the "Legendre transform" 

V($) =F(J($)) + $J($) 

and look at 

— — = + J + = J 

If I turn off the sources, this derivative vanishes. Thus the minimum of V tells us the 
expectation value of the field. This quantity V is the "effective potential." 
But now let me confuse you by looking at the second derivative of V 

d 2 V dJ 



a$2 gq, 
Actually, it is easier to look at the inverse 

d$ d 2 F 



dJ dJ 2 



(0 2 )-(0) 2 = ((0-(0)) 2 )>O. 



Thus this second derivative has a single sign! This shows we are actually looking for a 
minimum and not a maximum of V, but in addition it implies that V can only have ONE 
minimum! 

So what is going on? Are phase transitions impossible? The more you pull, the larger 
the expectation of $ should be. It won't go back. Physically, we must do Maxwell's 
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construction. If we force the expectation of to lie between two distinct stable phases, 
the system phase separates into a mixture of the two. Note that there is no large volume 
limit in the above discussion. However other definitions of V can allow a small barrier at 
finite volume due to surface tension effects. A mixed phase must contain interfaces, and 
their energy represents a barrier. 

First order transitions have a discontinuity in the internal energy, representing a latent 
heat. The barrier in the effective potential (modulo the above discussion) allows meta- 
stability and hysteresis. Water in a clean container can "bump" rather unpleasantly. 
Explosives last in a meta-stable state for long periods. Are things actually analytic as you 
pass through the transition? No, there is an essential singularity that I will now discuss. 
The free energies of the phases match at the transition; suppose I can expand as we go 
through it 

AF = C(p-p c ) 

Let the surface tension between the two phases be a. Creating a bubble of radius r costs 
free energy 

. 47rr 3 . ^ , 9 
E(r) = -— — AF + Anr 2 (j 

This has maximum energy at 

= Anr 2 AF + 8ivra 

or 

r = -2a I AF 
At this point the energy of the bubble is 

167ra 3 

' - max 



3(AF) 2 

As we approach the transition, this radius goes to infinity and this "semi-classical" argu- 
ment becomes rigorous. This energy represents a barrier to bubble nucleation, which is 
suppressed by the Boltzmann weight 

C 

P ~ exp(-(3E max ) = exp(- ^ ) 

An essential singularity appears in the physics as one passes into a meta-stable state. Even 
though things look very analytic, they are not. Since the meta-stable state can decay, this 
expression represents an "imaginary part" for the free energy of the unstable phase. 

Now for a brief discussion on some aspects of Goldstone Bosons. Suppose I have a 
conserved current 

d^jf, = 

so the corresponding charge Q = J d s xj (x) is a constant 

§ = -,[«,«] = a 
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Suppose, however, that the vacuum is not a singlet under this charge 

Q|0)^0 

Then there cannot be a mass gap in the theory. Consider the state 

exp(z# J d s xj (x)e- £x2 )\0) 

where e is a convenient cutoff and 6 some parameter. As epsilon goes to zero this state by 
assumption is not the vacuum, but the expectation value of the Hamiltonian goes to zero 
(normalize so the ground state energy is zero). Thus "spontaneously broken symmetries" 
have no mass gap, i.e. the theory contains states of arbitrarily low energy. These are 
manifested as massless particles called Goldstone bosons. 

Free massless field theory is a marvelous example of all this where everything can be 
worked out. The massless equation of motion 

d^d^cj) = 

can be written in the form 

d^U = 

where 

3 ix = 

The broken symmetry is the invariance of the Lagrangian L = J d A x{d lx (j)) 2 under shifts of 
the field 

(j) — > (j) + c 

Note that jo = do4> = 7r, the conjugate variable to <fi. One can work out explicitly 

(O|exp(z0 J d 3 x Jo (x)e-^ 2 / 2 )\0). 

but we can save ourselves the work using dimensional analysis. The field <p has dimensions 
of inverse length, while jo goes as inverse length squared. Thus 6 above has units of inverse 
length. These are the same dimensions as e 2 . Now for a free theory by Wick's theorem 
the answer must be Gaussian in 6>, so we conclude that the above overlap must go as 

exp(-C6> 2 /e 4 ) 

where C is some non-vanishing dimensionless number. This expression rapidly goes to 
zero as epsilon becomes small, showing that the vacuum is indeed not invariant under the 
symmetry. As e goes to zero, we obtain a new vacuum that is not even in the same Hilbert 
space. Its overlap with any polynomial of fields on the original vacuum vanishes. 
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Two dimensions give some interesting twists on this argument. Now the scalar field 
is dimensionless, and the current has dimensions of inverse length. Thus theta is dimen- 
sionless and we expect 

(O|exp(z0 J dxj (x)e- ex2/2 )\0) ~ exp(-C6> 2 ) 

There can be no log(e) since there is nothing to set the scale. Thus the vacuum is not 
invariant, but the symmetry relation does not give you a fully independent state. This is 
clearly a borderline case, and for an interacting theory the massless particles can be lost. 

This is tied in with the propagator in two dimensions not being a distribution. Put 
in a small mass cutoff. Then 



A(x) = 



d 2 p 



-ipx 



(2tt) 2 p 2 + m 2 



gets very singular as m goes to zero. Consider the simple test function e~ x I 2 . This is 
infinitely differentiable and well behaved at infinity. Now integrate this test function with 
the free particle propagator 



/ 



d 2 xA(x)e x2/2 ~ / I V ' „ ~ log(l/m) =^ m ^o oo 

J (27r) z p z + m z 



Thus A(x) is not a tempered distribution, contrary to the basic assumptions of quantum 
field theory. However Green's functions of j are tempered distributions since they involve 
derivatives that kill the divergent part. In most cases Goldstone bosons are lost in two 
dimensions, however, if they are free, as in the above case, they can exist. The X-Y model, 
with spins in U(l), has a massless phase, but no long range order. Lore is that higher 
symmetries only have massive phases, but Seiler and Patrascioiu have argued that this 
may be wrong. 

I now turn to the renormalization group, which I approach via the Migdal-Kadanoff 
approximate recursion relations. I start with a discussion of decimation. Let us go back 
to the Ising model in one dimension 

/ f3+H -/3 \ N 

Z = Tt {e-P e?- H ) =TtT 
Let me sum over every other spin, giving 

Z = Tr(T') N/2 

where 



r 



,2(fi+H) + e -2(3 e H + e -H 

e H +e -H e 2(fi-H) +e -2/3 
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We now match this with the original form of T 



/ e (3'+H' e -/3' \ 
T ' = C [ e -P' e P'-H> ) 

We see that exactly the same physics occurs on a lattice of twice the spacing and new 
couplings (/?', H'). The values of C, (3' and H' are fixed by the three equations 

Ce-P' =e H + e~ H 

Ce f3'+H' =e 2 { p+H) +e -2f3 
Ce^'- H ' =e 2(/3-H) +e -2/3 



This process is called decimation, i.e. integrating out some of the degrees of freedom. To 
simplify the equations, I turn off H, obtaining 

/?' = ^log(cosh(2/?)) 

This can be written in a form reminiscent of our earlier recursions 



The only fixed point occurs at (3 = 0. The new coupling is always less than the old one as 
long as beta is positive. Repeating this as an iteration drives any (3 to zero. 

It is instructive to extend this to non integer decimations. For this, write the transfer 
matrix in the form T ss / = cosh(/3)(l + ss' tanh(,3)). The above decimation by a factor of 
two involves the sum 



I ^(1 + s lS2 t)(l + s 2 s 3 t) = (1 + Sl s 3 t 2 ) 

S2 



or simply tanh(/3) — > tanh 2 (/3). Interpolate this to rescaling by a factor of 1 + A, taking 
tanh(,9) — > tanh 1+A (/?). Infinitesimally, this reduces to 

^ ^ — ~ a^- = cosh(,9) sinh(/3) log(tanh(/3)) 

This is the renormalization group equation for this system. The right hand side is negative 
for all positive f3. As the lattice spacing varies from zero to infinity, the coupling f3 flows 
from the ultraviolet fixed point at infinity to the infrared fixed point at zero. 
Suppose a system has a non-trivial fixed point satisfying 

af a =AGJ-ft). 
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Figure 9. Moving bonds around. 



This has the solution 

j3 = I3 C + Ca x 

Since a ~ l/£, this says 

l/£~(/?-/3 c ) lA 

This is the renormalization group way of seeing how non-trivial exponents can arise as one 
approaches a critical point. 

Going on to more dimensions we loose the exactness and must make approximations. 
Integrating out a site in more than one dimension introduces couplings between all sites 
to which it is coupled. Integrating the sites along a line couples all spins attached to that 
line. Integrating out all but the corners on a block requires an infinite number of couplings. 
This makes things less than rigorous, but can imagine a similar coupling "flow" in a higher 
space. 

Making an approximation by moving bonds around allows one to analytically study 
these flows. Imagine integrating out every other site in say the x direction. To avoid 
long range couplings being generated in the y direction, follow Kadanoff and "move" the 
y bonds to sites not being integrated over, as sketched in Fig. 9. Every second y bond 
becomes twice as strong, and then the earlier x decimation can be carried out on the 
remaining sites. Thus we relate the model at [3 x ,f3 y to that at 

/3; = ^log(cosh(2/3 x )) 

Now repeat this for the y direction. The resulting transformation is asymmetric due to the 
approximations. To get more symmetric, do things differentially, using the earlier equation 
for the x coupling and = f3 for the bond moving. The total change of coupling is then 

dd 

o-f- = cosh(/3) sinh(/3) log(tanh(/3)) + (d- 1)(3 
da 

Here I insert a factor of d — 1 to allow for bond moving in all directions but the decimation 
one. The result is exact in one dimension, and for d = 2 it still gives the exact f3 c . I plot 
this function in Fig. 10 for d = 2. 

The renormalization group relates theories with different lattice spacings. If we could 
keep track of an infinite number of couplings, the procedure would be "exact," but in 
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Figure 10. The Migdal-Kadanoff recursion relation for d = 2. 
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Figure 11. A generic renormalization group flow. 

reality we usually need some truncation. Continuing to integrate out degrees of freedom, 
the couplings flow and might reach some "fixed point." With two couplings, there can be 
an attractive "sheet" towards which couplings flow, and then they go towards the fixed 
point, as sketched in Fig. 11. If the fixed point has only one attractive direction, then two 
different models that flow towards that same fixed point will have the same physics. This 
is universality, i.e. exponents are the same for all these models with the same attractor. 

So mean field theory describes a phase transition in terms of a changing classical 
ground state as parameters are varied, and the renormalization group description is in 
terms of a flow through a complex coupling constant space. When should we rely on which 
picture? 



cosh(x)*sinh(x)*log(tanh(x))+x - 
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Some hints come from dimensional analysis, although, in ignoring non-perturbative 
effects that might occur at strong coupling, the following arguments are not rigorous. 

d — 2 

In d dimensions a conventional scalar field has dimensions of . Thus the coupling 

constant A in an interaction of form f d d x X(/) n has dimensions of M d ~ n ^~ . On a lattice 
of spacing a, the natural unit of dimension M is the inverse lattice spacing. Thus without 
any special tuning, the renormalized coupling at some fixed physical scale would naturally 
run as A ~ a n ^~~ d . As long as the exponent in this expression is positive, i.e. 



we expect the coupling to become "irrelevant" in the continuum limit. The fixed point is 
driven towards zero in the corresponding direction. If d exceeds four, this is the case for 
all interactions. (I ignore 4> 3 in 6 dimensions because of stability problems.) This suggests 
that four dimensions is a critical case, with mean field theory giving the right qualitative 
critical behavior for all larger dimensions. In four dimensions we have several possible 
"renormalizable" couplings which are dimensionless, suggesting logarithmic corrections to 
the simple dimensional arguments. Indeed, four-dimensional non-abelian gauge theories 
should display exactly such a logarithmic flow; this is asymptotic freedom. 

This simple dimensional argument applied to the mass term suggests it would flow 
towards infinity in all dimensions. For a conventional phase transition, something must 
be tuned to a critical point. In statistical mechanics this is the temperature. In field 
theory language we usually remap this onto a tuning of the mass term, saying that the 
transition occurs as some scalar mass goes through zero. This tuning of scalar mass terms 
required for a continuum limit seems unnatural and is one of the unsatisfying features of 
the standard model, driving particle physicists to try to unravel how the Higg's mechanism 
really works. 

Recently there has been considerable interest in statistical systems that become critical 
without any tuning of parameters. This phenomenon of "self-organized criticality" may 
explain the fractal structure of much of the world around us. However field theoretical 
applications of this concept remain elusive. 

In non-Abelian gauge theories with massless fermions, chiral symmetry protects the 
mass from renormalization, avoiding any special tuning. Indeed, these models exhibit the 
amazing phenomenon of dimensional transmutation: all dimensionless parameters in the 
continuum limit are completely determined by the basic structure of the initial Lagrangian, 
without any continuous parameters to tune. In the limit of vanishing pion mass, the rho 
to nucleon mass ratio should be determined from first principles; it is the goal of lattice 
gauge theory to calculate just such numbers. 

As we go below four dimensions, this dimensional argument suggests that several 
couplings can become "relevant," requiring the renormalization group picture of flow to- 
wards a non-trivial fixed point. Above two dimensions the finite number of renormalizable 
couplings corresponds to the renormalization group argument for a finite number of "uni- 
versality classes," corresponding to different basic symmetries. 

One might imagine dimensionality as being a continuously variable parameter. Then 
just below four dimensions a renormalizable coupling becomes "super-renormalizable" and 
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a new non-trivial fixed point breaks away from vanishing coupling. Near four dimensions 
this point is at small coupling, forming the basis for an expansion in 4 — d. This has become 
a major industry, making remarkably accurate predictions for critical exponents in three 
dimensions. 

Now I return to lattice gauge theory and discuss how pure glue at finite temperature 
mimics a three state model and might be expected to have a first order deconfining tran- 
sition. I lead into this with a bit of group theory. Consider some compact group G with 
elements g. There exists a unique measure 

J dgf(g) = J dgf{g x g) = J dgf{gg x ) 

where I normalize J dg = 1. (Non-compact groups might have different normalizations for 

left and right.) For example, with U(l) I can take g = e %e and then J dg = d9/2n. For 
SU{2), write g = cio + ia-a and then f dg = f d 4 a 5(d^ + a 2 — 1), i.e. the surface of a four 
dimensional sphere, an S3. 

This group integration extracts the "singlet" part of a function in the following sense. 
Suppose / is a "class function," i.e. f(g) = figiggi 1 )- Then I can expand it in traces 
over the various irreducible representations R of the group 

f(g) = ^2fRXR(g) 

R 

where the character xr{q) = T^^Rig) an d Mn(g) is the matrix representing g in represen- 
tation R. These representations include the trivial one, R = 0, the fundamental one R = F, 
the adjoint R = A, and generally infinitely many more. For irreducible representations, 
the characters x satisfy an orthogonality condition 

J d g x* R (g) XR'(g) = $rr> 

From this and Xo(g) = 1, we see 

J dgf(g) = f . 

If we insert a character in some other representation we obtain 

j dgx* R {g) f{g) = fR. 

This allows us to do some specific integrals integrals quite easily. For example, with SU(3) 
we have 

/ dg (Trg) 3 = 1 

JSU(3) 

since there is only one singlet in the famous decomposition 3®3(g>3 = 1©8©8©10. This 
integral lies at the base of the argument below for a first order chiral transition with three 
massless quark flavors. 
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So lets apply this to lattice gauge theory. On each bond of our hyper-cubic lattice we 
have a group element Uij. The Wilson action multiplies these around elementary squares 
and constructs Up = UxU^U^U^. The partition function is 



Formally, as argued earlier, this represents something like Z = Tre aHNt with a the temporal 
lattice spacing. The picture, however, is a bit more complicated due to gauge invariance. 
If we put a group element gi on each site, we can imagine taking Uj — > giUijgJ 1 . This 
change cancels from the action. 

Gauge invariance leads to the possibility of gauge fixing. This can be done much more 
generally, but for now suppose I forget to integrate over one link and define 



On the integrand I can do a gauge transformation and then use the invariance of measure 
to find 



So Z{Uq) doesn't depend on Uq. Since my measure is normalized, Z = Z . I can continue 
this and forget to integrate over more links. As long as no closed loops are fixed, then 
the partition function is unchanged. A closed loop is a gauge invariant observable, so we 
better not be able to fix it. 

In the temporal gauge, all time-like links are set to unity. In this gauge the above 
transfer matrix argument reduces the path integral to Hamiltonian lattice gauge theory. 
But on finite periodic temporal lattice, temporal links at a given spatial site form a closed 
loop. Thus one cannot gauge fix all of them. At each spatial location we must leave 
one temporal link unfixed, take it to be at time 0. What does integrating over this link 
correspond to? 

In Hamiltonian language, there is an operator Ri(g) that does a local gauge rotation 
at site i. In particular, for a link to a spatial neighbor 



These are all operators in the Hilbert space of the Hamiltonian approach. What the path 
integral formally reduces to is 



From the above discussion we see that this last integration projects out a gauge singlet. 
This operator imposes the lattice-gauge-theory version of Gauss's law 





p 



Z(U Q ) = Zigzag' 1 ) 



R i (g)U ij R- 1 (g)=gU 
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Now we can generalize and consider not projecting out the singlet everywhere. In 
particular, at one site I might want to put down a quark-like source. To do this I simply 
insert the character for the desired representation 



The ratio of the new partition function to the old is the Wilson line or the Polyakov loop. 
Going back to the path integral, it is just the expectation of a product of temporal links 
wrapping around the time direction. This Hamiltonian argument explicitly shows how it 
represents the energy carried by a fixed source in the fundamental representation of the 
gauge group. 

For the quark-less theory, the unfixed temporal links at time have a global symmetry 
under the center of the gauge group. For the SU(3) of the strong interactions the center is 
the set {1, e ±27ri//3 }. By definition, center elements commute with all group elements, and 
the global change gi — > e 27 ™/ 3 ^ will cancel out of the temporal plaquettes, each of which 
involves one gi and one g~ l . This is exactly the same symmetry as for the 3 state Potts 
model, which I argued above should have a first order phase transition. This prediction is 
well verified by numerical simulation. 

The quark kinetic term explicitly breaks this symmetry, so the transition might go 
away. For massless quarks, it is instead the global chiral symmetry that becomes relevant. 
This suggests second order for two flavors. For three light flavors the suggestion is first 
order since there is a quadratic term in the mean field expansion. This arises since the 
product of three fundamental representations contains a singlet piece, as we well know 
because three quarks can combine to form a gauge singlet baryon. Generalizing our earlier 
mean field equation to an SU (3) spin system gives 



where I use the earlier SU(3) example integral. The quadratic term means that the solu- 
tion jumps disco ntinuously, just as argued earlier for the Potts model. The interpolation 
between the small and large mass limits for various numbers of flavors is a major area of 
current study. 




3M = 



/ dg e 2dpMReTr 9Trg ^ d/3M + ( rf/ #M) 2 /2 + . . . 
J dg e 2d/3MReTr 3 ~ 1 + (d(3M) 2 . . . 



d/3M + (d/3M) 2 /2... 
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